JOURNAL OF APPLIED PHYSICS 


VOLUME 90, NUMBER 5 


1 SEPTEMBER 2001 


Impact of gas heating in inductively coupled plasmas 

D. B. Hash, 3 ' D. Bose, b) M. V. V. S. Rao, b) B. A. Cruden, b ' M. Meyyappan, 
and S. P, Sharma 

NASA Ames Research Center, Moffett Field, California 94035 
(Received 29 March 2001; accepted for publication 14 June 2001) 

Recently it has been recognized that the neutral gas in inductively coupled plasma reactors heats up 
significantly during processing. The resulting gas density variations across the reactor affect reaction 
rates, radical densities, plasma characteristics, and uniformity within the reactor. A self-consistent 
model that couples the plasma generation and transport to the gas flow and heating has been 
developed and used to study CF 4 discharges. A Langmuir probe has been used to measure radial 
profiles of electron density and temperature. The model predictions agree well with the experimental 
results. As a result of these comparisons along with the poorer performance of the model without the 
gas-plasma coupling, the importance of gas heating in plasma processing has been verified. 
© 2001 American Institute of Physics. [DOI; 10.1063/1.1390503] 


I. INTRODUCTION 

Inductively coupled plasma (ICP) sources have recently 
become popular in semiconductor processing. They provide 
a high plasma density ( 10 l °- 10 ,2 cm -3 ) and allow indepen- 
dent control of the plasma density and ion energy at the 
wafer surface. Plasma etching in an ICP reactor is a complex 
process characterized by strong interactions among plasma 
chemistry, surface chemistry, power coupling, plasma phys- 
ics, and gas flow. Plasma diagnostics as well as reactor mod- 
eling are important in developing an understanding of reac- 
tive discharges employed in plasma processing. Such an 
understanding is valuable in equipment and process design, 
and process control. 

Because of the low pressures at which ICPs and other 
plasma sources operate, simulation tools employed to ac- 
quire such insights may use a kinetic, fluid, or hybrid ap- 
proach. The first self-consistent two-dimensional nonequilib- 
rium simulations of a rf glow discharge were published by 
Wu employing both a fully fluid model 1 and a hybrid model 2 
with a Monte Carlo technique for electrons and a fluid rep- 
resentation for ions. Other particle electron-fluid ion simu- 
lations soon followed. The work of Ventzek et al . 3 laid the 
groundwork for the widely used Hybrid Plasma Equipment 
Model code. 4-10 Similarly, contributions were made to hy- 
brid research by Fiala el al. n During the same time period, 
Graves et al. ~ 14 introduced a particle ion-fluid electron hy- 
brid approach for simulating an electron cyclotron resonance 
plasma. A large number of research teams however have fo- 
cused on the development of a fully fluid approach. 1 15 33 

In the early 1990s when two-dimensional simulations of 
radio frequency plasmas began to appear, most of the afore- 
mentioned research, whether hybrid or fluid, neglected the 
coupling between the neutral gas and the plasma which re- 
cent experiments 34 have shown to be important. There were a 
few exceptions. Ventzek et al . 5 solved a decoupled hydrody- 
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namics model to provide the gas velocity flow field. Mass, 
momentum, and energy conservation equations were em- 
ployed but the energy equation did not include heating from 
inelastic collisions. Moreover, the only information used 
from the neutral gas simulation was the velocity field. They 
concluded that the effects of the gas flow on the global 
plasma parameters were secondary for low pressure systems. 
The Modular Plasma Reactor Simulator code of Lym- 
beropoulos and Economou 22 recognized the importance of 
gas-plasma coupling but only included mass conservation 
for the neutrals, neglecting momentum and energy conserva- 
tion. The convective term V ■ pu in the mass conservation 
was approximated by a relaxation term A p/r Tes . Finally, 
Fiala et al. 11 included mass and momentum equations for the 
neutral species but did not solve the neutral energy conser- 
vation. 

All of this research neglected what is now seen as a 
significant effect in plasma processing, that of neutral gas 
heating. Gas temperatures on the order of 1000 K have been 
observed in an argon ICP at 240 W, 35 1400 K in a 
C1 2 /BC1 j/N 2 plasma at 1000 W, 36 and 2100 K in a point-of- 
use plasma abatement tool at 1200 W. 37 This effect has only 
recently been addressed in the modeling community by 
Bukowski eta/., 38 Nam eta/., 39 and the cfdplasma code 
under the direction of Kolobov at CFDRC in Huntsville, and 
Kiehlbauch and Graves. 37 Bukowski et al. allow for gas 
heating through ion-neutral elastic and inelastic collisions, 
Franck-Condon heating, and ion-ion recombination. Nam 
et al. do include neutral gas conservation equations for mass, 
momentum, and energy. However, they make the assumption 
that the energy transfer from the plasma to the neutral gas 
through collisions is 5% of the deposited power in an ad hoc 
manner. Kiehlbauch and Graves present the most thorough 
representation to date of the true physics by modeling 
Franck-Condon heating, electron impact vibrational excita- 
tion of neutrals, gas phase exothermic reactions, ion-neutral 
gas heating through inelastic and elastic collisions, and ion- 
ion recombination. They show that in fact the power trans- 
ferred to the neutrals in their CF 4 work is more on the order 
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of 25%-30% of the rf power. But this of course is dependent 
on the gas itself and the reaction mechanisms involved. Re- 
gardless, it is now evident that detailed modeling of the ac- 
tual energy transfer mechanism is necessary for accurate 
simulations. 

The purpose of the present research is to further the un- 
derstanding of the importance of gas heating in plasma pro- 
cessing and also further validate the Semiconductor Equip- 
ment Modeling Software (sems) code under development at 
NASA Ames Research Center. CF 4 mixtures have been ex- 
tensively investigated in the Gaseous Electronics Conference 
(GEC) reference cell but none have been published in the 
ICP mode. The following sections describe the computa- 
tional model and the experimental measurements of an 
ICP-CF 4 plasma in a GEC reference cell. Comparisons of 
the sems predictions to Langmuir probe measurements of 
electron temperature and density will then be used as the 
metric for code validation and for the the discussion of the 
importance of gas heating. Further experimental results and 
analysis are presented in two companion papers. The first 40 
characterizes the electron energy distribution function, and 
the second 41 uses Fourier transform infrared spectroscopy to 
characterize the degree of CF 4 dissociation which is therein 
compared to results from this work. 


II. COMPUTATIONAL MODELING 

sems consists of a comprehensive fluid model for simu- 
lating high-density plasma reactors. The model couples the 
plasma and neutral gas transport and heat transfer in a self- 
consistent manner. No assumptions regarding the energy 
transferred from the plasma to the gas are made. Instead, 
similar to Kiehlbauch and Graves , 37 the energy exchange 
mechanisms between the plasma and the gas are modeled. 
Details regarding the model have been presented 
previously 42 44 but will be briefly addressed here. The model 
solves a coupled set of equations including mass conserva- 
tion for each species, a total momentum equation, and ther- 
mal energy conservation equations for the neutrals and elec- 
trons. 

A. Species continuity 

^ + V P ,u=-V J, + 2 4. (1) 

0 1 r 


where the subscript i represents a single ion species, the 
number density n„ is for a single neutral species in m 3 , the 
polarizability a n is also for a single neutral species with units 
of A 3 , and the reduced mass m r is measured in g/mole. The 
ion mobility is then determined from summing over all ion- 
neutral interactions with Blanc's law 



(3) 


where X n is the mole fraction of the neutral species n . 

The subscript r in the source term of the continuity 
equation denotes an individual reaction with the term R SI 
representing the production of species s from reaction r de- 
fined as follows: 


R^MsiK-Okrll ‘ r , 

k r =AT a e\p( — E a /T). (4) 

M is the molecular weight, and k r is the reaction rate con- 
stant in an Arrhenius form where E a is the activation energy 
and T is the effective temperature for the reaction, i.e., ion, 
electron, or neutral temperature depending on the reaction. 
All reactions used in the simulation are fit to this form. When 
available, reaction cross sections are integrated over an as- 
sumed Maxwellian distribution to provide the reaction rates. 
These reaction rates are then curve fit to the Arrhenius form. 
The terms v' and v" in the above equation are the stoichio- 
metric coefficients for the reaction 


2 2 v 'Lxs , (5) 

s s 

where \ is the chemical symbol. 


B. Total momentum 

dpM 

+ Vpuu= - Vp+ V • t. (6) 


Here, p is the total pressure, and r is the shear stress. 


C. Neutral thermal energy 

+ u=-V- 2 (q J + M*)-/>nV-u 

01 s , neutrals 


Here, the subscript ,v denotes an individual species and 
p s its density, u is the mass averaged velocity of the mixture, 
and J is the diffusion mass flux, p s \ s , with V, representing 
the species diffusion velocity. The diffusion fluxes are deter- 
mined from the self-consistent effective binary diffusion for- 
mulation of Ramshaw and Chang . 45 47 For the ions, this re- 
quires knowledge of the mobilities p, to determine the 
diffusivities D through the Einstein relation, p=\e\ !kT,D, 
where e is the elementary charge. The mobilities are deter- 
mined from the Langevin formula 48 

3.73X 10 22 , 

P,, n n n = — (nr/Vs), (2) 

\ a n m r 


+ 2 ~ VPt + (2elastic+ (2 ce 

5 , neutrals f's 

+ 2 2 RsAK ( 7 ) 

r J t neutrals 

Here the subscript n denotes a property of all the neutrals as 
a combined set, e n is the neutral thermal energy, e n 
= 2 S ±i,ePsCv,sT, q is the heat flux, and h is the specific 
enthalpy. The terms V • p„e n u, V q, and V • h J are the trans- 
port of energy due to convection, thermal conduction, and 
diffusion, respectively. The next two terms in the equation 
are due to pressure volume work. The remaining terms rep- 
resent sources to the neutral thermal energy equation 
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whereby the heating of the gas is modeled. The G elastlc term 
models the energy transfer due to elastic collisions between 
electrons and neutrals. The Q cc term models the energy trans- 
fer due to charge exchange collisions between ions and neu- 
trals. The R ST AH r term models the chemical energy released 
or absorbed in a reaction where A H° r is the heat of reaction. 
This term then includes gas heating from ion-ion recombi- 
nation, dissociative recombination, gas phase exothermic re- 
actions, and associative detachment. 

Viscous dissipation, neutral vibrational excitation, and 
the Franck-Condon effect are not considered. Kiehlbauch 
and Graves 37 show viscous dissipation to be negligible, sems 
has the capability of modeling vibrational excitation but this 
was shown only to represent 10% of the heating of the neu- 
tral gas 37 and is therefore also ignored. The Franck-Condon 
effect is a result of the quantized energy states allowable to a 
molecule. When an electron impacts a molecule because of 
the Franck-Condon rule there are only certain states possible 
for the molecule to achieve. If dissociation occurs as a result 
of the electron impact, this can result in extra energy being 
deposited into the molecule fragments resulting from the dif- 
ference in the energy required for dissociation and the al- 
lowed excited state. Information concerning the additional 
energy funneled to the fragments is largely unavailable and 
thus has not been included in the present work. 

D. Electron thermal energy 

do e e e i e 

— ~ + V-p e e e u= -V ■(<\, + hJ e )-p t! V-u+ — -Vp e 

(it f) e 

(2 elastic , s ^er A// r 

.v * e r 

+ 2)o e b- ( 8 ) 

Here, the term jjgEg is the power deposited into the plasma 
from the coil assuming fully collisional heating, j e is the 
azimuthal plasma current, and Eg is the azimuthal electric 
field. 

Maxwell’s equations are solved separately to compute 
E g and j e where j e =aE ti using Ohm’s law with u being the 
plasma conductivity. The Maxwell’s equations are then 
coupled at regular intervals to the conservation equations. 
For these results, a Poisson solver is not used, and the space 
charge induced electric fields are given by the Boltzmann 
relation assuming that the plasma is quasineutral. The most 
recent version of sems 49 includes Poisson’s equation and 
separate ion momentum equations to consider the important 
inertial effects at low pressure, but these capabilities are not 
utilized in the present work. 


III. cf 4 modeling 

Fluorocarbon plasmas such as CF 4 , CHF 3 , and C 4 F 8 are 
used in oxide etching. Given the significance of this process 
to the microelectronics industry, there has been an emphasis 
on basic research to understand the characteristics of fluoro- 
carbon plasmas. One of the main impediments to modeling 


has been either the lack of measurements of reaction cross 
sections or uncertainty in those that have been measured. For 
CF 4 , this is especially relevant for the reactions involving 
negative ions. 50 The reaction set used in this work is summa- 
rized in Table I, consisting of electron impact, neutral- 
neutral, and neutral-ion exchange reactions. The rates are 
given in Arrhenius form and are determined from integration 
over a Maxwellian distribution of the reaction cross sections 
when available. Construction of the reaction set in Table I 
assumes a negligible concentration of negative ions in a CF 4 
plasma and avoids the uncertainty in the electronegative re- 
action mechanisms. Experimental confirmation for the large 
presence or absence of negative ions in CF 4 plasmas is 
scarce. In our work, a mass spectrometer with a positively 
biased extraction orifice was mounted on the side of the 
chamber to detect negative ions. While this approach dis- 
turbs the plasma and may not be representative of the true 
negative ion density, the negative ion flux was 2 orders of 
magnitude lower than the positive ion flux, suggesting the 
negative ion density is fairly small. For completeness. Table 
II lists additional reactions which need to be considered if F 
is found in significant amounts; these reactions include dis- 
sociative attachment, polar dissociation, associative detach- 
ment, and ion-ion recombination. The effect of the inclusion 
of this reaction set is discussed in the results. 

All the simulations are for a power input of 3(H) W with 
an estimated power absorption of 80% based on current- 
voltage (/-V0 probe and reflected power measurements. The 
CF 4 flow rate is 20 seem, and the wall temperatures are fixed 
at 400 K. The ion temperature is specified at 0.2 eV for use 
in calculating the ion diffusivities. Three different pressures 
are considered; 10, 30, and 50 mTorr. The fluorocarbon gas is 
modeled with 13 species: CF 4 , CF-„ CF 2 , CF, F, C, F 2 , CF,’ , 
CF 2 , CF + , F + , F“, and e. All ions are neutralized upon 
impact with the wall, and atomic fluorine is allowed to re- 
combine on the walls with a probability of 0.02. 51 The trans- 
port properties for the neutrals are determined using the 
Lennard-Jones parameters shown in Table III and the tech- 
niques in Gupta et al. 52 The ion transport properties are de- 
termined using the polarizabilities given in the same table, 
and the electron transport properties are given from integrat- 
ing the electron momentum transfer cross section presented 
in Christophorou et al . 53 


IV. EXPERIMENT 

Our experimental setup consists of a GEC reference cell 
modified for inductive coupling shown in Figs. 1-5. A de- 
tailed description of the apparatus, the Langmuir probe (LP) 
construction, and measurement procedures has been recently 
presented 54 and hence, only a brief account is provided here. 
The GEC cell is powered by a flat five-turn and 3.5 in. diam. 
inductive coil made of a 0.125 in. diam copper tube. The 
inductive coil is mounted at the top of the chamber and is 
concentric with the center of the chamber to maintain dis- 
charge symmetry. The coil is insulated from the plasma by a 
quartz window of 0.375 in. thickness and 5 in. diameter. The 
water-cooled lower electrode (4 in. diameter) is made of 304 
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TABLE I. CF 4 reaction mechanism. Species: CF ( , CF 3 , CF 2 , CF, C, F, F 2 , Cl ,' , CF 2 , CF ' , F + , e. 


Reaction 

Rate coefficient (cm 3 s ‘) 

A//” (eV)‘ 

Reference 

e 1 CF 3 -.CF 3 1 2e 

1.213X 10“ 9 7 99076 exp{— 9.2727/7",) 

9.09 

Christophorou et alP 

e + CF 2 — >CF 2 + 2e 

4.071 X 10“ 9 7 01029 exp(— 9.4004/7,) 

11.5 

Christophorou et a l. 53 

e + CF— >CF + + 2e 

3.62 1 X 1 O' 9 7° 8665 exp( - 1 1 . 1 98/7",) 

9.17 

Christophorou et a/. 53 

e + F— »F 4 + 2e 

4.702 X 10~ 9 I*''* exp( - 17.586 IT,) 

17.5 

Meyyappan anti Govindan 50 

e + CF 4 — >CF 3 +F+ 2e 

4.153X 10' 9 7" 33554 exp(- 15.599/7",) 

14.8 

Christophorou el a/. 53 

f + CF„~>CF 2 +2F+2e 

2. 1 04 X 10- 9 7 9 * 103 exp( -23. 164/7",) 

20.8 

Christophorou et a/. 5 ’ 

e + CF 4 — *CF + + 3F+ 2e 

2.729 X 10~ 9 7"“ 7747 exp< -33. 101/7,) 

23.9 

Christophorou et al. 5> 

e + CF 3 — >CF 2 +F+2e 

9.001 X 10“ 9 T? 4111 exp( - 17.463/7",) 

15.2 

Christophorou et a/. 53 

e + CF 3 — *CF + +2F+ 2<? 

5.045 X 10^ 9 7 05737 e xp( -21.174/T,) 

18.2 

Christophorou et alP 

e + CF 2 — >CF + +F+ 2e 

4.737 X 10“ 9 7, 06696 exp( - 14.3 10/7",) 

14.6 

Christophorou et al* 3 

e + CF 4 — *CF 3 +F+* > 

3.257X 10“ 9 7 06900 exp( - 14.370/7",) 

5.67 

Christophorou et a/. 53 

<? + CF 4 — *CF 2 +2F+ e 

2.221 X 10- 10 T °; 9939 exp(- 14.775/7",) 

9.32 

Christophorou et alP 

e -f CF 4 — »CF+3F+e 

7.289X lO" 9 r;° 1653 exp( -28.261/7",) 

14.7 

Christophorou et alP 

e + CF 3 — >CF 2 +F + e 

3.257 X 1 0 " 9 7° 6906 exp( - 3.65/7",) 

3.65 

Christophorou et al . 53b 

e 4- CF 3 — >CF+2F+ € 

2.221 X 10“ 10 7? 9939 exp(- 9.04/7",) 

9.04 

Christophorou et al S3b 

e + CFj— *CF+F+e 

3.257 X 10" 9 7° 6900 exp( -5.39/7",) 

5.39 

Christophorou el al.* 3 * 

e + CF— >C+F + f 

3.257 X 10“ 9 T 96906 exp( -5.69/7,) 

5.69 

Christophorou et al. 5Sb 

e -t F 2 — »F+F+ e 

3.931 X 10~ 9 7 32270 exp( -2.544/7,) 

1.65 

Edelson and Flamm 56 

e + CF 3 ->CF 2 + F 

6X 10 * 

-5.44 

Nighan and Wiegand 57c 

e + CFj — >CF+F 

6X 10~ s 

-6.10 

Nighan and Wiegand S7c 

e + CF + — <C+F 

6X10"* 

-3.47 

Nighan and Wiegand 5 " 9 

cf,+f 2 -*cf 4 +f 

3.882X 10- 12 # ex p(-0.1256/7) 

-4.03 

Modica ajid Sillers 58 

CF 2 +F 2 — CFj+F 

3.542X10-“ \jf exp(- 0.092/7) 

-2.0 

Modica and Sillers 5 " 

F + +CFj--.F+CFj 

1X10 -9 

-8.4 

Xu et alP 

f + +cf 4 -*f 2 +cf 3 + 

1X1Q- 9 

-4.37 

Xu et alP 


*Heat of reaction. 

b Estimate based on CF 4 analogy. 

‘Estimate based on similar systems. 


stainless steel and is grounded in the present study. The in- 
ductive coil is connected to a 13.56 MHz rf power supply 
through a matching network. 

The LP measurements were made using a commercially 
available rf compensated Smart Probe made by Scientific 


Systems Inc. The body of the probe is made of a ceramic 
tube with a tapered ceramic tip holder. A platinum wire 380 
/xm in diameter and 2.5 mm in length is chosen for the probe 
tip. The probe is mounted on a motion feedthrough and in- 
serted into the plasma through one of the radial ports. The 


TABLE II. Additional reactions involving negative ions, F . 


Reaction 

Rate coefficient (cm 3 s *) 

A H\ (eV)‘ 

Reference 

e + CF 4 — *CF 3 +F 

1.603X 10~ 10 7," 04289 exp( -4.139/7,) 

2.21 

Christophorou et al. 53 

e + CF 3 — >CF 2 +F 

1.603 X 10" 10 7-; 0 - 42119 exp( -0.18/7,) 

0.18 

Christophorou et al 53b 

e + CF 2 -*CF+F~ 

1.603X 10- 10 r; 042 " 9 exp( - 1.93/7,) 

1.93 

Christophorou et al.* 3b 

e + CF— >C4 F" 

1.603X 10" 10 r; 042 * 9 exp( -2.22/7,) 

2.22 

Christophorou et a!. 5it 

e -t F 2 — >F+F _ 

3.974 X 10- 9 T~ 33759 exp{ -0.2256/7,) 

- 1.82 

Hayashi and Nimura 60 

e + CF 4 -♦CF, 3 -t-F" + e 

5.263 X 10 - 12 7? 7677 exp( - 14.687/7,) 

11.3 

Christophorou et al . 

c + CFj— >CF 2 +F" + c 

5.263 X 10- 12 T? 7677 exp( - 1 1.678/7,) 

11.7 

Christophorou et al. 53b 

e + F” — »F+ 2e 

5.456X 10- 10 7 3l75 exp(- 3.296/7,) 

3.47 

Edelson and Hamm 50 

CF 4 +F - — >CF 4 +F+e 

2.132X 10-* 7 3131 exp(- 1.448/7,) 

3.47 

Peko el al . 01 

CF,+F-— CF 4 + e 

5X10- 10 

-2.21 

Nighan and Wiegand 57c 

CF 2 +F _ — »CF 3 +e 

5X10" to 

-0.18 

Nighan and Wiegand 57c 

CF+F~— »CF 2 + e 

5X10" 10 

- 1.93 

Nighan and Wiegand 57c 

C+F - — tCF+e 

5X10-'° 

-2.22 

Nighan and Wiegand 57 ‘ 

F+F~— >F 2 + e 

1.4X10- 10 exp(- 1.9/7,) 

1.82 

Mandl 62 

cf;+f-cf 3 +f 

8.8X 10" 7 

- 5.62 

Hebner mid Miller''' 

cf 2 * +f _ -.cf 2 +f 

8.8 X 10- 7 

-8.03 

Hebner and Miller 63 

CF + + F ->CF+F 

8.8X10" 7 

-5.70 

Hebner and Miller 63 

F 3 + F-^F + F 

4.OXIO- 7 

- 14.0 

Nighan and Wiegand 579 


■“Heat of reaction. 

b Estimate based on CF 4 analogy. 

‘Estimate based on similar systems. 
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TABLE HI. Parameters for transport properties. 


Lennard-Jones parameters 


Molecule 

Vanishing distance, 
,T (A) 

Well depth, 
elk (K) 

Polarizability*, 

a (A 5 ) 

cf 4 

4.662 

134 

2.82 

cf 3 

4.320 

121 

2.82 b 

cf 2 

3.977 

108 

2.82 b 

CF 

3.635 

94.2 

l.8 c 

F 

2.968 

112.6 

0.6 

C 

3.385 

30.6 

1.5 

F 2 

3.357 

112.6 

0.9' 1 


"See Ref. 48. 

b Estimate based on ( ' F, analogy. 

‘Estimate based on CO analogy. 

d Estimate based on relation between polarizability of other diatomics and 
their atoms. 

motion of the probe is provided by a computer-controlled 
stepper motor. The probe tip is cleaned between each data 
point collection by positively pulse biasing the probe tip so 
that the tip becomes white hot and desorbs deposited species. 
Data are taken at seven radial locations from the center of the 
chamber to 6 cm in 1 cm increments. The probe was axially 
positioned in the center of the plasma between the quartz 
window and the lower electrode. 

The effective electron temperature and density are deter- 
mined by the appropriate weighted integration of the electron 
energy distribution function (EEDF), calculated from the 
second derivative of the probe l-V characteristic. 54 The 
EEDF is discretized at 40-50 points, each of which is aver- 
aged over 1000 measured data points. The uncertainty of the 
LP measurements is estimated to be 15% for the electron 
density and ±0.5eV for the electron temperature. 

V. RESULTS AND DISCUSSION 

A. Plasma characteristics 

To illustrate the plasma characteristics, two-dimensional 
contour plots of the properties throughout the reactor for the 



FIG. 1. Contours of electron number density ( 10 1 ' cm ’) at 50 mTorr from 
sems simulation. 



lation. 

50 mTorr case are presented. Figure 1 displays the electron 
number density profiles in the ICP reactor at this pressure. 
The peak density is 4.25x 10 11 cirT 5 in the center of the 
plasma. An off-axis maximum is not predicted in the simu- 
lation, nor is it observed in the experiments (compare Fig. 7). 
Figure 2 shows the electron temperature. The toroidal shape 
of the temperature profile also seen in the experimental mea- 
surements corresponds to the toroidal power deposition of 
the ICP. 

Figure 3 illustrates the contours of the CF 4 mole fraction 
which show significant consumption of the feed gas in the 
plasma chamber through various electron impact reactions. 
Figure 4 displays the CF 3 density indicating higher concen- 
trations of the radical near the quartz window and the lower 
electrode. These higher concentrations result from electron 
impact dissociation of CF 4 and neutralization of the domi- 
nant ion, CF 3 , at the surfaces (see Fig. 5). CF, is formed in 
the chamber by electron impact of the feed gas as well as the 



FIG. 3. Contours of CF 4 mole fraction at 50 mTorr from sems simulation. 
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FIG. 4. Contours of CF 3 number density ( 10 13 cm 3 ) at 50 mTon from sems 
simulation. 


CF 3 radical and then neutralizes upon impact with the quartz 
window and lower electrode. The peak CF 3 density in the 
central region of the plasma corresponds to the peak electron 
density in Fig. 1. The predictions reveal that CF 3 is the 
dominant ion. 

B. Experimental comparison 

Figure 6 compares the sems results using the reaction set 
in Table I to the radial Langmuir probe measurements for 
electron temperature made halfway between the quartz win- 
dow and the lower electrode. The results are plotted from the 
axis of the reactor to a distance of 8 cm from the axis. The 
lower electrode extends out to 5.08 cm. Included with this 
comparison are drift-diffusion results where gas flow and 
heating are ignored. For the latter, the neutral gas thermal 
energy equation and the total momentum equation are re- 



FIG. 5. Contours of CF 3 number density ( 10 11 cm 3 ) at 50 mTorr from 
sems simulation. 




(b) 



FIG. 6. Comparison of electron temperature from Langmuir probe measure- 
ments (at the axial mid-plane) and simulations: (a) 50 mTorr, (b) 30 mTorr, 
and (c) 10 mTorr. Simulations using approximate drift-diffusion formulation 
(dashed lines) and the complete model including gas flow and heating (solid 
lines) are both compared against experiment (circles). 


moved, and the species continuity equations are replaced by 
the drift-diffusion assumption in which convection is ne- 
glected. In addition, the total density, gas pressure, and gas 
temperature are assumed constant throughout the reactor for 
the drift-diffusion case. It is seen from the figure that the 
electron temperature rises with a decrease in pressure as ex- 
pected. In all three cases, the full model including the effects 
of gas flow and gas heating predicts the experimental results 
better than the drift-diffusion approximation. Inclusion of the 
gas heating effect in the full model is primarily responsible 
for this predictive capability. Note that, unlike in the larger 
commercial reactors with substantial flow rates, the pressure 
drop across the GEC cell is minimal and hence its impact on 
gas flow and gas density variation is negligible. Therefore, 
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(a) 




FIG. 7. Comparison of electron numher density from Langmuir probe mea- 
surements (at the axial mid-plane) and simulations: (a) 50 rnTorr, (b) 30 
mTorr, and (c) 10 mTorr. Simulations using approximate drift-diffusion for- 
mulation (dashed lines) and the complete model including gas flow and 
heating (solid lines) are both compared against experiment (circles). 

the increased gas temperature due to charge exchange colli- 
sions and dissociative recombination directly results in a de- 
crease in neutral gas densities and in general in a variation of 
neutral gas density across the reactor. The decrease in gas 
density results in lower electron impact production rates, 
leading to a reduced electron density and also results in an 
increase in electron temperature, relative to the drift- 
diffusion approximation with its fixed gas temperature and 
neutral gas density assumptions. 

Figure 7 compares the SF.MS results to the Langmuir 
probe measurements for electron number density. Once again 
the drift-diffusion results are included. The predictions from 
the full model of sems compare reasonably well with the 
experiments and better than that of the drift-diffusion ap- 



proximation. The experimental results for the 10 mTorr case 
fall midway between the predictions of the full model and 
the drift-diffusion model. As the collision frequency de- 
creases with pressure, gas heating becomes less important at 
low pressures resulting in less difference between the full 
model and the drift-diffusion model for the 10 mTorr case 
than the higher pressure cases. The reduction in collision 
frequency means less exothermic reactions that funnel en- 
ergy from the plasma to the gas. The decrease in gas tem- 
perature with pressure is clearly seen in Fig. 8 which dis- 
plays the radial profiles of the gas temperature from the 
center of the plasma to the reactor wall along the axial plane 
located at the midpoint of the plasma. The peak neutral gas 
temperature falls from 550 K at 50 mTorr to 490 K at 30 
mTorr and then to 460 K at 10 mTorr. Also, in Fig. 7, note 
that the drift-diffusion model consistently predicts higher 
densities than the full model (and the experiment) which is 
an artifact resulting from the fixed lower neutral gas tem- 
perature (400 K) and the resulting higher neutral density in 
the drift-diffusion approximation as discussed earlier. 

Although in some low aspect ratio ICPs off-axis maxi- 
mums in electron density can be seen, this is not the case 
here in either the experimental measurements or the sems 
predictions. The simulations do a good job of predicting the 
trends of the density profiles including on-axis maximums. 
The experimental comparison of electron density with pre- 
dictions is not as good as for the electron temperature. This is 
often the case in modeling given the uncertainty in the power 
absorption efficiency. It is well known that electron tempera- 
ture does not vary much with power while electron density is 
strongly power dependent. The assumption of 80% power 
absorption we have used may not be completely accurate, 
resulting in the discrepancies. In fact, later measurements 
suggest that it is closer to 70%. Another possible reason is 
that the electron impact rate constants (in Table I) were ob- 
tained by integrating the cross sections using a Maxwellian 
EEDF while the probe measurements 40 indicate a non- 
Maxwellian EEDF. 

C. Effect of negative ions 

All of the simulations discussed above were also re- 
peated with the inclusion of negative ions (F~) and the cor- 
responding reaction set in Table II. Figure 9 shows the effect 
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FIG. 9. Effect of negative ions on electron temperature predictions at 50 
mTorr. 


of including the electronegative reactions on electron tem- 
perature for the 50 mTorr case. The sems results including 
the effects of gas flow and heating show an increase in elec- 
tron temperature, although the results still fall within the ex- 
perimental error. This increase is a result of the increased 
heating from the highly exothermic reactions of ion-ion re- 
combination, which efficiently pump energy from the plasma 
into the neutral gas. As described previously, elevated gas 
temperatures result in lower neutral densities, which result in 
lower electron impact production. This produces a lower 
electron density, and a higher electron temperature. The drift- 
diffusion model, with the fixed gas temperature approxima- 
tion, does not model this effect and thus only exhibits slight 
differences with the inclusion of negative ions; the small 
change seen is solely due to differences in radical densities 
resulting from the electronegative reactions. 

Figure 10 illustrates the effect of including the electrone- 
gative reactions on the neutral gas temperature. The gas tem- 
perature now reaches a maximum of 670 K compared to 550 
K for the case without negative ions. Moreover, the profiles 
are different. Without negative ions, the peak heating is out- 
side the active plasma region where significant dissociative 



FIG, 10. Effect of negative ions on gas temperature (K) at 50 mTorr. 


recombination takes place. In contrast, the peak heating is 
seen at the center of the chamber when negative ions are 
included due to the highly exothermic ion-ion recombina- 
tion reactions. The rate of these reactions in the central 
plasma volume is high due to the high local densities of both 
positive and negative ions. 

The 50 mTorr case predicts a plasma with the ratio of 
number of negative ions to electrons less than one ( n In,, 
< 1 ). However, the 30 mTorr case predicts a plasma with a 
ratio greater than one («_/«,,> 1) in the plasma chamber. 
This highly overpredicted negative ion density results in an 
erroneous electron temperature mid-plasma radial profile 
with an axis value of 10 eV, more than twice that of the 
experimentally measured value. The authors are lead to the 
conclusion that the electronegative reaction mechanisms are 
suspect. This was further evidenced when it was found that 
the 10 mTorr case was not possible to simulate with the sems 
code using the electronegative reaction set. Due to the over- 
prediction of negative ions, it was impossible to maintain the 
plasma. 

In the experimental work, biasing on the Langmuir 
probe could induce instabilities, suggesting the conditions 
are close to an unstable operating regime. Under some ex- 
perimental conditions (typically higher pressure and lower 
power), the discharge was in fact seen to extinguish or 
switch to capacitive coupling mode. This type of instability 
has been associated with the presence of negative ions, 55 and 
is not inconsistent with the effect of high negative ion con- 
centrations observed in the sf.ms code. However, the sems 
predictions in this regard are not completely consistent with 
the experimental data as it is possible to obtain operating 
conditions at 10 mTorr and 300 W without this instability, 
and the instability is greater at higher, not lower, pressures. 
This discrepancy in trends further calls into question the ac- 
curacy of our electronegative reaction set. 

The experimental measurements presented in Fig. 6 may 
indicate some presence of negative ions. Otherwise, it would 
be expected that the toroidal shape of the electron tempera- 
ture profile would dissipate with decreasing pressure as the 
electron thermal conductivity correspondingly increases. As 
a result, the off-axis maximums in the radial electron tem- 
perature profile should become less pronounced and move 
toward the axis. Instead, in the experiment, the off-axis 
maximums may be maintained because of localized produc- 
tion of negative ions due to the toroidal deposition as shown 
in Fig. 11. The resulting localized depletion of electrons 
could sustain the off-axis maximums in the electron tempera- 
ture. Further experimental and computational investigations 
into the presence of negative ions in pure CF 4 plasmas are 
warranted and will be the subject of future work. 

An alternate and/or complementary explanation may 
arise from the fact that the EEDF is non-Maxwellian con- 
trary to the assumptions of our model. It has been shown in 
research by Kortshagen et al. 6 * that the EEDF has some 
bearing on the radial profiles of electron temperature. This 
may also contribute to the discrepancies in the radial trends 
predicted by sems and measured experimentally. 
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FIG. 11. Contours of F number density ( 10" cm 7 at 50 mTorr from sems 
simulation. 

VI. CONCLUDING REMARKS 

The demands in equipment and process design, arising 
from the continuous scaling down of device feature sizes, 
scaling up of reactor dimensions, and from the increased 
costs of new fabs, have made modeling and diagnostics valu- 
able in integrated circuit manufacturing. Development of 
technology computer aided design tools require accurate 
models that capture the process physics and chemistry. In 
plasma processing, models have historically done well in 
capturing the discharge physics and power coupling features 
but overlooked issues related to the neutral gas. In the 
present work, we have used the sems code, which includes 
gas How and gas heating models coupled to plasma genera- 
tion and transport, to study a ICP-CF 4 discharge in a GEC 
reference plasma reactor. The analysis shows significant 
heating of the gas and temperature variations across the re- 
actor which affect the plasma characteristics and radical den- 
sities. We have also measured the electron density and tem- 
perature in the GEC reactor using a Langmuir probe. The 
predictions agree well with the measurements. When the full 
model in sems is relaxed to the familiar drift-diffusion for- 
mulation, by dropping the Navier-Stokes and gas-energy 
equations and assuming a constant gas density and tempera- 
ture throughout the reactor, the model results depart more 
from the measurements. Finally, the experiment could indi- 
cate some presence of negative ions and/or the need for more 
accurate modeling of the EEDF through the use of a Boltz- 
mann solver. This would likely improve the quality of our 
comparisons to the experiment in regard to the prediction of 
the off-axis maximums in electron temperature. However, 
current knowledge of negative ion reaction kinetics is sus- 
pect, and further detailed research into the actual electrone- 
gative reaction mechanisms and cross sections is warranted. 
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